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When are active Brownian particles and run-and-tumble particles 
equivalent? Consequences for motility-induced phase separation 
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Abstract -Active Brownian particles (ABPs, such as self-phoretic colloids) swim at fixed speed v 
along a body-axis u that rotates by slow angular diffusion. Run-and-tumble particles (RTPs, such 
as motile bacteria) swim with constant u until a random tumble event suddenly decorrelates the 
orientation. We show that when the motility parameters depend on density p but not on u, the 
coarse-grained fluctuating hydrodynamics of interacting ABPs and RTPs can be mapped onto 
each other and are thus strictly equivalent. In both steeply enough decreasing v(p) causes 

phase separation in dimensions d = 2,3, even when no attractive forces act between the particles. 
This points to a generic role for motility-induced phase separation in active matter. However, 
we show that the ABP/RTP equivalence does not automatically extend to the more general case 
of u-dependent motilities. 
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The physics of self-propelled colloidal particles repre- 
sents a central focus of research into 'active matter' [l}{3]. 
In active matter, a continuous supply of energy destroys 
microscopic time-reversibility and allows phenomena to 
arise that are impossible in thermal equilibrium systems, 
where detailed balance restores time-reversal symmetry in 
the steady state. 

One class of self propelled colloids is represented by 
motile bacteria such as Escherichia coli, which move in 
a sequence of 'runs' - periods of almost straight-line mo- 
tion at near-constant speed v - punctuated by sudden and 
rapid reassignments of direction, or 'tumbles', occurring 
at random with rate a [4j[5]. (Certain algae have sim- 
ilar motion [6j.) At large scales, and under conditions 
of constant v,a, a single run-and-tumble particle (RTP) 
performs a random- walk of diffusivity D Q = v 2 /ad (d is 
dimensionality) which cannot be distinguished from the 
equilibrium dynamics of non-interacting passive Brownian 
particles (PBPs). Nonetheless, it was shown previously 
that if their swim speed v(p) decreases fast enough with 
the local particle density p, RTPs undergo phase separa- 
tion between a dense, slow-swimming fluid and a dilute 
fast-swimming one [7j. This happens even in the com- 
plete absence of interparticle attractive forces, which are 
a prerequisite for fluid-fluid phase separation in thermal 



PBPs. Coupled to population dynamics, this motility- 
induced phase separation was implicated in the patterning 
of growing bacterial colonies [8 10. 



A second class of self-propelled colloids are represented 
by synthetic (bi-)metallised colloids which, usually by 
catalysing the breakdown of hydrogen peroxide, create a 
self-phoretic local chemical motor [IT- 15 . These differ 
from RTPs in that the direction of swimming only changes 
gradually, by rotational diffusion (with angular diffusivity 
D r ). This rotation is usually Brownian in origin, hence 
the name ABPs (active Brownian particles). Biological 
ABPs also exist, including non-tumbling E coli mutants 
and other organisms that self-propel without tumbling. 
Clearly an isolated ABP again performs a random walk 
(Do oc v 2 /D r ) at large scales. 

When all microscopic parameters are uniform and 
isotropic, the large scale dynamics of RTPs, ABPs and 
PBPs are thus described by the diffusion equation. Equat- 
ing their diffusivities thus trivially makes their dynamics 
identical. In the converse case, for instance when interac- 
tions, external potentials, or spatially varying swimming 
parameters are present, the dynamics become more com- 
plex and typically involve a drift term: 

p = -V • [V(v, a, D r ,...)p- D(v, a, D r ,... )Vp] (1) 
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While the microscopic parameters can clearly still be cho- 
sen so that the diffusivities of the three processes are 
equal, there is no reason why the equality of the drift 
terms V should then follow as a by-product. Following 
Schnitzer [2] , this point can be understood by considering 
PBPs in which either viscosity or temperature are non- 
uniform in space, both leading to the same nonuniform 
D(r). As required by detailed balance for an isothermal 
system, V is zero in the first case despite the nonuniform 
diffusivity. But if the same diffusivity variation arises in- 
stead through a non- uniform T(r), detailed balance is bro- 
ken, V is nonzero, and no generic Einstein-like relation re- 
lates D to V. Likewise, because ABPs and RTPs are out 
of equilibrium, the mere existence of a diffusive behaviour 
in the uniform case does not establish the equivalence of 
ABPs with either PBPs or RTPs in more general situa- 
tions. It is therefore intriguing that recent simulations of 
spherical ABPs with purely repulsive interactions showed 
clear signs of phase separation 16 17 , something very 



18 



similar to what is observed for RTPs 

By analogy with what we know for RTPs, such a phase 
separation is explicable if the repulsions act to effectively 
decrease v at high density (by caging, for instance) . How- 
ever, the growing interest in ABPs [3j (exemplified by sev- 
eral other recent theory /simulation studies and continuing 
experiments 19 -26]) clearly requires a treatment of this 
subject that lies beyond mere analogy. Are ABPs and 
RTPs actually equivalent, and if so under what conditions? 

This Letter provides a partial answer to that important 
question. At time scales large compared to the angular 
reorientation time r (~ 1/D r or ~ l/cn), and length scales 
large compared to I = vt, we establish a leading-order 
exact equivalence between ABPs and RTPs so long as the 
motility parameters v, a, D r are independent of particle 
orientation u. This result, which extends to the case 
where a Brownian translational diffusivity D t is added, 
also holds when these parameters depend (slowly) on po- 
sition. Following [7j, this allows treatment of interacting 
particles, as long as the effect of interactions can be faith- 
fully described by density-dependent swimming parame- 
ters. This can for instance be the case when interactions 
are mediated by rapidly diffusing chemicals [9 10 or when 
mean-field treatments are valid (either as a result of molec- 
ular chaos or when interactions are averaged over many 
neighbouring particles as in [18]). Thus the ABP-RTP 
equivalence indeed holds under the conditions (generalised 
here to d > 1) shown in 7 to cause motility-induced bulk 
phase separation in RTPs. Accordingly ABPs will gener- 
ically show the same sort of phase separation, subject to 
the same requirement of a sufficiently decreasing v(p). 

Because we assume slow variations of all quantities on 
the scale of £, we cannot rule out differences between 
ABPs and RTPs in situations involving ordering at shorter 
scales, such as the solidification transition addressed in 



23 . Moreover, as we prove later on, no general equiv- 
alence exists for cases involving u-dependent motility pa- 
rameters. (This shows that while some degree of simi- 



larity between ABPs and RTPs might be expected, any 
actual equivalence between them is not automatic.) In- 
stances of such anisotropic dynamics include bacterial rec- 
24p7 , as modelled by an orientation-dependent 



tification 
tumble rate 

speed acquires a gravitational component 14||22l|28 and 



28 , sedimentation, where the particle 



systems of ABPs in which interactions promote align- 
ment [3{ [2T|[29}|33] . Furthermore, hydrodynamic interac- 
tions, neglected entirely here, have subtle consequences 
[34p8] that might well differ between RTPs and ABPs. 
Despite these caveats, it is interesting to establish an 
asymptotic equivalence between ABPs and RTPs in a sec- 
tor of parameter space that includes a nontrivial many- 
body phenomenon, namely motility-induced bulk phase 
separation. 

Fluctuating hydrodynamics. We start by analysing 
a single particle undergoing a generalized run-and-tumble 
dynamics with both rotational and translational diffu- 
sions, and later recover RTP and ABP as limiting cases. 
We allow, v, a, D r and D t to depend on position but not 
orientation. The probability density ^(r, u, t) of finding a 
particle at position r moving in direction u obeys exactly 
(in d = 2, 3) 

y)(r,u) = -V.[vuV>(r,u)] + V u [L> r V u ^(r, u)] 

a f (2) 
+ V(Z>tW(r, u)) - aip(r, u) + - / ^(r, u')dtf 

where V u is the rotational gradient acting on u and the 
integral is over the unit sphere |u'| = 1 of area J7. The first 
term on the right is the divergence of the advective current 
resulting from self propulsion and the last two terms are 
loss and gain terms due to tumbling in and out of the 
direction u. We next decompose ip as 

ip(r, u, t) = ip + p.u + Q : (uu - l/d) + Q[ip] (3) 

Here, ip, p, Q are functions of (r, t) but not u, which pa- 
rameterize the zeroth, first and second angular (g? = 2) 
or spherical (d = 3) harmonic components of ip, while 6 
projects onto the higher harmonic components [39] . 
Integrating eq. |2]) over u gives 



V> = -^V(up)+V(AVp) 



(4) 



Using the orthogonality of angular/spherical harmonics, 
one can then obtain equations order by order. Multiplying 
eq. |2]) by u and integrating over u then yields, after some 
algebra |40], for each component of p 



Pa = -V„(«¥») ~ (D r (d - 1) + a)p a + V{D t Vp a ) 
— B abcd \7 b (vQ cd ) 



(5) 



where B abcd = (S ac S bd + S ad S bc ~ 2S ab S cd /d)/ (d + 2). 

Similarly, multiplying eq. ([2| by uu — l/d and integrat- 
ing over u yields (with \abc defined below) 

Qab = - < ^T^B abcd V c (vp d ) - V c (Xabc) ,„s 

2 (6) 
- (2dD r + a)Q ab + V(AVQ a6 ) 
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Equation Q states that the rate of change of the prob- 
ability density of finding the particle at position r irre- 
spectively of its orientation (zeroth harmonic, tp) is the 
negative divergence of a flux J = vp/d— D t Vtp which de- 
pends solely on the two first harmonics. In addition to 
the diffusive flux — _D ( Vy>, there is a contribution of the 
self-propulsion along p. 

In eq. ^ , one sees that p relaxes via a flux term but also 
directly by angular diffusion and tumbling: here (1 — d) is 
the eigenvalue of A u for the first harmonic (likewise — 2d 
for the second in eq. ([6| ) . The last term gives the contribu- 
tion to p arising from the Q term within V.[uin/;] in eq. ([2]). 
There are no contributions to p from the higher harmon- 
ics but there is one from the zeroth harmonic: the V a (v<p) 
term. It encodes the fact that, although an isotropic dis- 
tribution of swimmers cannot transport density if D t = 0, 
it can still create anisotropy; for instance if the swim speed 
or density is higher to the right of the origin, an initially 
isotropic density at the origin soon develops an excess of 
left-moving swimmers. 

A corresponding set of remarks apply to eq. ^ for the 
time evolution of Q. Here B aoc jV c (vpd) is a flux contri- 
bution from the first harmonic whereas Xabc (whose form 
we don't need) arises from the higher harmonics in Q[ip]. 
Finally, the evolution of all remaining harmonics is ob- 
tained, if needed, by projecting eq. Q using 0. 

So far, beyond the assumed isotropy of v(r), D r; t(r) and 
a(r), no approximation has been made; eqs. (|4]|6| are ex- 
act results for the time evolution of the zeroth, first and 
second harmonics of tp(r,u,t). At this stage, despite 
similar structures, the dynamics of ABP and RTP cannot 
be mapped onto each other. Equating the prefactors of 
p a in ([5]) indeed requires (d — l)D r = a which makes the 
prefactors of Q a b in ([6| unequal since then 2dD r 7^ a. 

We now wish to coarse-grain these equations to obtain 
a 'diffusion/drift' expression for the flux J that involves 
only the conserved probability density <p(r,t), the slowly 
varying parameters D rt ,v,a, and their first derivatives: 



ip = -V.J ; J = V(r,t)<p- D(r,t)V(p 



(7) 



Following |7j, we use a gradient expansion and first note 
that ip is the only slow mode: its relaxation time is of 
order ~ (V) -1 whereas all other harmonics relax in times 
of order ~ 1. We thus assume p = Q = ©[^] =0 when 
calculating the quasi-stationary current J. By itself, this 
creates a non-local constitutive relation between J and 
ip which still involves all harmonics. Next we carry out 
explicitly the gradient expansion, yielding for Q: 



" " ~ 2dD r + a 



Qij 



C(V 2 ) 



from which the quasi-stationary p then follows via Q as 
1 



(d-l)D r + a 



V(t^) + 0(V 2 ) (8) 



Fortunately therefore, to diffusion-drift order, closure is 
achieved without needing further information on harmon- 
ics beyond the first. To this order eq. ^ holds, with 



d(d- l)D r + da 



V(vip) - D t Vtp 



so that the diffusivity and drift velocity obey 



D = 



d(d - l)D r + da 



-D t ; V = 



-vVv 



d(d - l)D r + da 



(9) 



(10) 



Eqs. (JtJ) and ( 10 1 give the evolution for the probability 
density of one particle at diffusion-drift level. They are 
equivalent to an Ito-Langevin equation for an individual 
particle position r /t (t) . Following [7] we can therefore now 
consider an assembly of particles whose motility parame- 
ters v, a, D r and D t depend on position through a set of 
(smooth) density functionals. The coarse-grained density 
p(r, t) on which these depend then obeys the many body 
Langevin equation [7] 



p = -V. (V[p]p - D[p]Vp + (2Dpf' 2 A 



(11) 



with white noise (A^r, t)Aj(r', i')) = %<5(r - r')8(t - t'). 
This noise term serves as a reminder that p(r, i) is not 
the one-particle probability (p(r,t) but a smooth coarse- 
graining of the collective density field J^fj. $( r ~ r M W) (with 
p a particle index) which evolves stochastically. The func- 
tionals v[p], a[p] and D t<r [p] in ( |To| then defines for the 
interacting particle system the many-body drift velocity 



and diffusivity V[p] and D[p] for use in eq. (11). 

We can now use eqs. ( 10|11 ) to compare the dynamics of 
RTPs and ABPs. From the expression of D and V in ( 10 ), 



we see that the tumble rate a and rotational diffusivity 
D t enter only through the combination [d — l)D r + a = 
a. Hence the pure RTP and pure ABP limits, alongside 
anything in between, are made equivalent by a suitable 
choice of a[p]. The effect of nonzero D t is the same in 
each case, although we note that, in many experimental 
cases, the self-propelling speed v is large enough that the 
translational diffusivity D t is negligible. (For instance in 
wild-type run-and-tumble bacteria, v 2 /(da) is about three 
orders of magnitude larger than D t .) In this limit, pure 
ABPs (D t = a = 0) and pure RTPs (D t = D r = 0) are 
equivalent up to the mapping (d — l)D r o a. 

Motility- induced phase separation. The many- 
body physics predicted by eqs. ( |10|11 ) for the large class 
of active particle systems addressed above essentially co- 
incides with that studied in [7J for RTPs in d — 1. (Note 
that ABPs cannot be defined in d — 1.) As shown there, 



because ( 11 ) is equivalent to a set of interacting PBPs with 



the specified diffusivity and drift functionals, a condition 
can be found under which the interacting active system 
will behave just like a fluid of PBPs with free energy (in 
thermal units) J-[p] = J~ e x[p\ + /p( m P ~ l)dx. The re- 
quired condition is [7] 

V([p],r)/D([p],r) = ~W(ST ex [p]/5p(r)) (12) 
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where the right hand side represents the force (i.e., excess 
chemical potential gradient) on a particle at r. 

In the case where D t = 0, the left hand side of fi"2| ) 
is simply — Vlnw[p] and we then require SF^p]/ 5p(r) = 
\n(v([p];r). When the swim speed depends only locally 
on density, so that w([p];r) = v(p(r)), we have T ex = 
/ fex(p(r))dr where f ex = J Q P In v(s)ds. (Note that when 
adD t is a nonzero constant, this result generalizes to 
fex = \ J P ln[ w ( s ) 2 + adD t ]ds.) This form of excess free 
energy leads at mean-field level to a spinodal instability 
whenever dv/dp < —v/p, and to binodal conditions for 
phase coexistence given by the usual common tangent con- 
struction on / = f ex + p(hip — 1) [7]. Beyond mean field, 
the steady state probability for density fluctuations is gov- 
erned by the 'Boltzmann-like' distribution exp[— F] and 
these fluctuations will eventually cause phase separation 
at all densities between the binodals. 

To understand the phase-separation dynamics in detail, 
to address the interfacial tension between phases, and also 
to confirm that the phase separation is equilibrium-like 
(rather than having, say, an ever-moving interface between 
the phases), would require a detailed examination of gra- 
dient terms lying beyond the present study. The same 
was initially true of the analysis of RTPs made in [7] , but 
the relevant gradient terms were later identified explicitly 
18 . However, some of the gradient terms neglected via 
eq. ^ could in principle differ between ABPs and RTPs, 
even if these share the same functional V[p] and D[p\. 
Accordingly, the somewhat technical exploration of these 
higher gradient terms is, for the ABP case, deferred to 
future work. 

The similarity between RTPs and ABPs undergoing 
phase separation can nevertheless be explored numerically 
by comparing simulations of large populations of interact- 
ing ABPs and RTPs. To do so we set D t = and con- 
sider a swimming speed w[p(x)] = vq exp[— \<p arctan(p/</>)] 
which decreases exponentially at low density before satu- 
rating at finite but non-zero swimming speecQ for larger 
p. Numerically, p(r) is computed by convoluting the 
number density 5Z M #(r ~ r /i) with the function /(r) = 
Z exp[— 1/(1 — r 2 /w 2 )] with Z a normalisation constant. 
For such a u[p(a;)], the free energy is bimodal as soon as 
\<j) > 2. Starting from identical uniform initial condi- 
tion, the spinodal decomposition of 6400 RTPs and ABPs 
is shown in figure [l] and is consistent with the aforemen- 
tioned mapping between ABPs and RTPs. 

We can now turn back to the simulations of 



16 17 



which show that excluded volume interactions promote 
phase separation of ABPs. This result, like that found 
numerically for RTPs on a lattice 18 , fits well with 



the scenario presented here, even though the interactions 
in 16 17 are collisional rather than via a smooth density 



functional. Notably, the collisions considered in 16 17 



do not directly affect the orientation u; at high density a 




, , , , 
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Fig. 1: Spinodal decomposition of 6400 ABP (left) and RTP 
(right) in 2D with v = 10, a = D r = 1.25, A = 0.45, </> = 7, 
w = 3.5 at times (top), 100 (middle) and 1000 (bottom). 
Each particle appears as an arrow pointing along its body axis 
and with a colorcode equal to the local density p(r). 



particle undergoes many collisions before one-body rota- 
tional diffusion relaxes its swimming direction. Modelling 
collisions via v(p) seems reasonable in this case. Indeed, 
v(p) was reported in 16 to decay linearly, a case shown 
in [18] to cause phase separation for RTPs. Although a 



microscopic model starting directly with an effective v(p) 
is only a mean-field-like description of the ABP simula- 
tions of 16,17 , the mechanistic similarity between their 



1 As in [7] , the density in the dense phase diverges in the absence 
of excluded volume interactions if v(p) vanishes at large p. 



phase separation and that of RTPs is now clear. 

Anisotropic motilities. We now turn to more general 
cases where the dynamics is not isotropic and discuss the 
limits of the equivalence between ABPs and RTPs. There 
are many ways of breaking isotropy and, as we show below, 
the mapping between RTPs and ABPs can be violated 
qualitatively, quantitatively, or preserved, depending on 
details of the microscopic dynamics. 

One simple way to remove isotropy is to add a uni- 
form field (such as gravity) that superimposes a Stokes 
drift velocity w on the self-propulsion. The overall ve- 
locity of a particle is then vu + w. This only affects the 
advection term in |2| and adds — V(w<p), —V b(wbPa) and 
— V c {w c Qab) to the right hand side of equations (4|6 1. The 
large time/scale expansion then yields Q ~ V as before 
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and one gets for p and ip 



1 



(D r (d-1) 



ip = -- viup) - v(-wip) ^ 

a 



V(D t Vv>) 

+ C(V 2 ) 

(13) 

At this coarse-grained level, there is again an exact equiv- 
alence between RTPs and ABPs upon identifying a and 
(d— l)D r . This is particularly interesting since a uniform 
field can induce polar order 22 , which thus does not by 



itself breaks the mapping between ABPs and RTPs. 

Another way of breaking isotropy is to have the swim- 
ming parameters depend on the local orientation u. For 
simplicity, we will only consider a first order spherical 
harmonic dependence D r (u) = Dq + • u and a(u) = 
ao(u) +ai-u. We will also take D t = for sake of clarity. 
Equation ^ now becomes 

ijj = -V.[vui/>] + V u • [(D^ + T>1 ■ u)V„V] 

1 (14) 
— (a + oli ■ u)ijj + a p + ,ol\ ■ p 

d 

Once again, we project this equation onto the three first 
harmonics. Tedious but straightforward algebra, that will 
be detailed elsewhere, show that equation Q is not mod- 
ified whereas ( 5|6 1 become 



Vi(vtp) - [a + (d- l)Dl\pi - ax,i<p 
B i3 ki[Vj{vQ k i) + aijQki + dD{ :Q k i\ 



(15) 



and 



d + 2 



Bijktf 'k[vpi] - ^ J tilS\f k S\^ n Vk{vTi mn ] 
2a + MDo 



d + 2 



-Qij - Bijke(ai,k + dD x ^)p t 



K> + 2(d + i) Dl , lB „„/d<>s, 



(3) c (3) 
Iran 



where = UiUjU k - (ui5 jk + Uj8 ik + u k 5ij)/{d + 2) 
is a tensor of 3rd order spherical harmonics and Tij k the 
component of ip along this tensor. 



Equation (15) already shows that whatever the choice 



of coefficients, there can never be an equivalence between 
ABPs and RTPs in this context, since the tumbling yields 
a contribution — ai t i<p, which has no counterpart stem- 
ming from the rotational diffusion and which does not dis- 
appear in the large time/scale limit. Interestingly, the ad- 
dition of a well chose angular potential for the rotational 
diffusion partially solves this problem. Such a potential 
produces an angular drift term ip ~ — V u (V u V>)j choos- 
ing a drift V u = — V u -D r then effectively replaces the 
term V u (£>(u)V u V>) in JO) by A u (D(u)ip), which adds a 
— (d — l)Di t iP term to (15). If this is done, the hydrody- 



namic equations for ABPs and RTPs do have similar struc- 
tures, but it is still not possible to choose motility param- 
eters to make all the coefficients match. Thus the map- 
ping breaks down quantitatively. While these technical 



subtleties will be presented in more detail elsewhere, the 
results above are enough to prove that no general equiv- 
alence between ABPs and RTPs exists, once anisotropic 
dynamics is present. 

Discussion. We comment finally on the role of the 



effective free energy ( 12 ) in relation to the concept of 
effective temperature. Several works have shown that 
active baths can be described by effective temperatures 
when dealing with passive tracers 41-43 or external po- 
tentials 14 28 , which implies that active particles may 
be considered as 'hot colloids'. For interacting ABPs, 
however, the usefulness of such effective temperatures has 
been questioned [16j[23]. We have shown above that the 
phenomenology expected for interacting ABPs is funda- 
mentally different from that of PBPs at any temperature. 
Nevertheless, under the specific conditions addressed here 



for which (12) is obeyed, equilibrium-like concepts do ap- 



ply: the steady states of ABPs then satisfy detailed bal- 
ance and so admit a free energy description. Currently we 
do not know the conditions (if any) under which a path 
can be traced from the microscopic dynamics of interact- 
ing ABPs directly to an effective free energy, avoiding the 
coarse graining steps that led us from ^ to (12). More 
generally, the validity or otherwise of effective temperature 
descriptions may largely depend on whether detailed bal- 
ance emerges in a coarse grained description even though 
it is, by definition, absent microscopically in all active sys- 
tems. 

Conclusion. In summary, this Letter establishes that, 
under conditions where their local swim speeds and reori- 
entation rates do not depend on the swimming direction, a 
large class of active particle models obey a common equa- 
tion at the diffusion-drift level. The ABP and RTP limits 
are connected by a mapping, a f-> (d — l)D r , between 
tumble rate a and rotational diffusivity D r . Were tum- 
bles to be considered of finite duration, one would expect 
a more complex mapping in which the swim speed v of the 
equivalent ABPs also depends on a [7|. 

The diffusion-drift description treats all nonconserved 
density harmonics as quasi-stationary, and expands the 
total particle flux to first order in spatial gradients; it 
therefore requires slow variations in the conserved density 
(or zeroth harmonic). Since a description at this level is 
sufficient to explain the motility-induced phase separation 
of RTPs whose swim speed is density-dependent, ABPs 
will show such phase separations under exactly the same 
conditions on v(p) as derived previously for RTPs. (In 
passing we have clarified that those conditions apply in 
d = 2,3 as well as for d = 1 as derived originally [7].) 
Motility-induced phase separation is caused by a feedback 
between (a) the tendency for particles to move slowly in re- 
gions of high density and (b) the tendency for particles to 
accumulate where they move slowly. The second of these 
is prohibited by detailed balance in thermal equilibrium. 
Our work shows that such phase separation (derived pre- 
viously only for RTPs in ID) is a robustly generic feature 
of interacting active particles. 
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Furthermore, in cases where the dynamics is not 
isotropic, we have shown that the mapping can some- 
times survive (even in the presence of polar order) but can 
also break down because of new couplings between har- 
monics of different orders. Nonetheless, progress might be 
possible when the angular dependencies on u are in turn 
expandable in low order angular or spherical harmonics as 
was the case in 



32 



Although the approach presented in this letter already 
can treat the sedimentation 14 22 and rectification 27 



problems alluded to previously, direct interparticle forces 
remain to be explored. For example, interparticle attrac- 
|25l|44] or nematic interactions 21 33 promote clus- 



tions |25[|44| or nematic interactions 
tering tendencies that seem distinct from the bulk phase 
separation caused by motility modulations alone. In gen- 
eral, interparticle forces create local orientational variation 
of the particle speed (just as the uniform force of sedi- 
mentation creates a global variation) [7] which effectively 
v(u) or v = v(p). Accordingly, despite some 



means v 



initial steps that have been made recently [22p9p0p2 



an 



explicit treatment of that case appears necessary as a first 
step towards a general dynamic density functional theory 
for interacting ABPs. Such a theory, which might involve 
analysis connected with that done for liquid crystals in 



45 , is a worthy goal for future research. 
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